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1. Introduction 

To relate observables measured in lattice QCD to their physical counterpart in the continuum, 
renormalisation is needed. Besides non-perturbative renormalisation also perturbative approaches 
are useful. In addition, it is useful to know as precisely as possible perturbative contributions to 
lattice observables assumed to show confinement properties in order to separate non-perturbative 
effects (condensates etc.). The gluon and the ghost propagator belong to these observables. 

It is well known that lattice perturbation theory (LPT) is much more involved compared to 
its continuum QCD counterpart. The complexity of diagrammatic approaches increases rapidly 
beyond the one-loop approximation. By now only a limited number of results up to two-loop 
accuracy have been obtained. 

Applying the standard Langevin dynamics [JT], |2[] to the problem of weak coupling expansions 
for lattice QCD, a powerful numerical approach for higher loop calculations - called numerical 
stochastic perturbation theory (NSPT) - has been proposed in Amongst other results un- 
quenched Nf = 2 Wilson loops up to 3-loop order plaquettes up to 8-loop order in pure 
QCD ^ and renormalisation constants related to the QCD pressure [^] have been calculated. 
There is ongoing progress in calculating high-loop perturbative renormalisation constants 0. As a 
new application we report here on a higher-loop calculation of the perturbative contributions to the 
gluon propagator in Landau gauge. For a similar related study see [Jsj] . More detailed results will 
be presented elsewhere. 

2. The Langevin equation 

The basis of stochastic quantisation is the Langevin equation derived from the Euclidean action 
that generates a (quasi-continuous) sequence of Euclidean field configurations. For gauge theories 
some special aspects related to the gauge redundancy have to be taken into account. Let t the 
Langevin time, then the Langevin equation reads 

^U x ^{t-r,) = i [V X ^S G [U] - r] x ^t)} U^farj) , (2.1) 
where -q = ^ a rj a T a is a random noise field with a Gaussian distribution satisfying 

= , (vt,(t)v b y ,v(.t'))^ = 2S ab V 5 xy 5(t - . (2.2) 

The notation (. . .) denotes an average over the (external) Gaussian stochastic measure. As for all 
Gaussian processes, higher cumulants vanish. T a are the (anti-hermitian) generators of the gauge 
group SU(N). The differential operator V XjAt = ^ a T a V" is the left Lie derivative for any 
function on the group and a partial derivative with respect to the links of the lattice. 

It can be proven that the gauge fields, in the limit of large t for the continuous-time Langevin 
equation, are distributed according to the Gibbs measure P[U] oc exp(— Sg[U]). In practice, 
the Langevin equation is solved by discretisation of time, t = ne, with running step number n. 
Therefore, in order to extract correct physical information, it is not only necessary to go to large t, 
but also to do the extrapolation e — » 0. For the solution of the Langevin equation we use the Euler 
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scheme in a way that guarantees all the links U x „ E SU(N) not to leave the group manifold: 



U X)fl {n + 1; rj) = exp(iF^ [C/, 77] ) U XtfJi {n; rj) 
F x ,n[U, rj\ = e V Xifl S G [U] + \/e • 



(2.3) 
(2.4) 



For £g[^] we take here the standard plaquette Wilson gauge action. 

For the stochastic perturbation theory it is substantial to consider each link matrix as an ex- 
pansion in the bare coupling constant g. Since = 2N/g 2 , the expansion reads 



1 + 



£ 

i>i 



(3- l/2 U«Uf,v). 



(2.5) 



It simplifies matters if one rescales the time step e = 0e. Upon the expansion in g, the Langevin 
equation transforms (2.3) into a system of simultaneous updates that takes the following form in 
terms of the expansion coefficients of U X:fl 



(2.6) 



and of the force F X)(t ( [2.4] ) 
U (1) (n+l) = U {1) {n) - F {l \n) 

U (2) (n + l) = U i2 \n)-F^(n) + ^(F {1 \n)) 2 -F^(n)U W {n) 



The random noise r/ enters only the lowest order equation through F^\ the lowest part of the force 



(2.4) analogous to the expansion (2.5). Higher orders are stochastic only by the noise fed in from 



the lower order terms. 



A similar expansion like (2.5) exists also for the (anti-hermitian) vector potential living in the 
algebra su(N), 

A x+fl/2 Jt; V )^Y, P~ l,2A %,2^ V) ■ ( 2 - 7 ) 
l>l 

Since the vector potential A x+ ^/2,ii is related to the links U xfJi via Aj. + £/ 2)/L4 = \ogU x ^, the sep- 
arate orders AW can be expressed via the orders XJ^ k \ Enforcing unitarity of the originally unex- 
panded U XifJi link fields is tantamount to enforcing anti-hermiticity and tracelessness of all orders 
of A x+ ^/ 2}tl . Whenever we speak about contributions of some order to an observable this has to be 
understood in the sense of an expansion 

(O) -» ^/T* /2 (0 (/) ). (2.8) 

l>0 

3. The (perturbative) gluon propagator 

The lattice gluon propagator D?~,(q) is the Fourier transform of the gluon two-point function, 
i.e. the expectation value 

DfM = (A;(k)A b u (-k)) = 5 ab D^ v (q) , (3.1) 

which is required to be color-diagonal. Here A^(k) is the Fourier transform of A" + ^/ 2 ^, and q 
denotes the physical discrete momentum 

g M (fc M ) = - a sin (j^j = 2 - sin (^) , k^ e (-L M /2, L„/2] (3.2) 



3 



The lattice gluon propagator in stochastic perturbation theory 



Arwed Schiller 



corresponding to the integer- valued Fourier momentum 4-vector k on the finite lattice. Some values 
q 2 (k) of the lattice momentum squared (the lattice equivalent of q 2 in the continuum limit) can be 
realized by different integer 4-tuples fo, &3, ^4). 

Assuming reality of the color components of the vector potential and rotational invariance of 
the two-point function, the continuum gluon propagator has the following general tensor structure 

D^{q) = U M „ - JD{q) + ^- , (3.3) 

with D(q 2 ) and F(q 2 ) being the transverse and longitudinal propagator, respectively. The longitu- 
dinal propagator F{q 2 ) vanishes in the Landau gauge. 

The lattice gluon propagator D^(g) depends on the lattice four-momentum q. Inspired by the 



continuum form (3.3) we consider the following lattice scalars ^ v q^D^(q)q v and ^ D^(q) 
that should survive the continuum limit. The first scalar vanishes exactly in lattice Landau gauge. 
In this gauge the second scalar function, corresponding to the tranverse part of the gluon propagator 



in the continuum limit, is denoted by 

i 

3 



1 4 

D{q) = -Y J D^{q). (3.4) 

On the lattice, this function shows the lower symmetry of the hypercubic group in that it depends 
on the scalar quantities X)u^ n > n = 1? 2, . . . rather than being a smooth function of only q 2 . 
Multiplying with q 2 or (aq) 2 we get the two versions of dressing functions or form factors 

Z(q) = q 2 D(q) , Z(aq) = (aq) 2 D(q) . (3.5) 



Using the expansion (2.1) and - /2 =J2 a T a A a ^~ /2 we obtain the different loop orders 



n of the perturbative gluon propagator (restricted to even powers of I in the sense of (2.8)), 



2n+l 



S^D^iq) = (J2 [^ {i) {k)A b ^ 2n+2 ^(-k)\ \ . (3.6) 

Note that already the tree-level contribution, Dj$, arises from the quantum fluctuations of the 
gauge fields with i = 1. Therefore, the tree-level result for the dressing function, Z^°\q) = 1 in 
the limit e —> for all sets of lattice momenta (k±, k%, ks, k^), is nontrivial and is obtained as the 
result of averaging. 

4. Practical implementation of NSPT 



Solving the coupled system of equations ( |2.6| ), one generates a configuration sequence of 



expanded gauge fields at finite e which can be used to measure the perturbatively constructed ob- 
servables. To study the limit e — > 0, we used e = 0.07, 0.05, 0.03, 0.02, 0.01. It is expected that the 
autocorrelation time r for a chosen observable extending over subsequent configurations increases 
with decreasing e. As reasonable compromise between computer time and autocorrelation we have 
measured the gluon propagator after each 20th Langevin step. The remaining autocorrelations are 
taking into account in the error estimate. 
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To obtain infinite volume perturbative loop results at vanishing lattice spacing, different lattice 
sizes have to be studied in addition. We have used L = 6,8, 10, 12(16) and studied the maximal 
loop number for the propagator n max = 4 (1). After reaching equilibrium for the largest Wilson 
loops to all orders, up to 60000 Langevin steps have been used to obtain up to 3000 measured 
gluon propagators. We checked that expectation values for odd powers of I for all observables 
indeed vanish in the ensemble average. 

The Landau gauge for all considered orders of the perturbative gauge fields (/ m ax = 10) is 
defined by the condition 



(0 

x+fi/2,^ 



A 



(0 

x—(i/2, fi 



(4.1) 



For the configurations used in measurements it is reached by an iterative gauge transformation 
using the expanded variant of Ui ifJ , = G{x) U x ^G\x + fi) . The gauge transformation is chosen 
as a perturbative variant of Fourier acceleration [P] with an optimal a = l/p 2 



f 

'max 



exp 



"2 
p2 



(4.2) 



p 2 is the non-zero eigenvalues of lattice — d 2 and F (F ) denotes the forward (backward) Fourier 
transform. The iterative procedure to reach the Landau gauge stops when, for all orders I, 



within double precision. 



In the course of the Langevin process we mainly follow the prescription given in pj]: after 
each Langevin step we perform a gauge transformation with (x) = exp — e d^A$ ilx and 
subtract zero momentum gauge field modes to all orders. This keeps the gauge field components 
finite. 



5. Selected results for the gluon propagator 



The measured gluon propagator data have been averaged over equivalent 4-tuples of lattice 
momenta and linearly extrapolated to the limit e = 0. In Figs. [T] we present the extrapolation at 
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Figure 1: Tree level (left) and one-loop (right) dressing function Z^ 0,1 ) (q) vs. q 2 at L = 16. 

lattice size 16 4 for the tree and one-loop dressing function Z^ 0,1 ) as function of q 2 defined in (3.2), 
together with original data at e = 0.07 and 0.03 . 
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Figure 2: One-loop dressing function Z^\aq) vs. (aq) 2 at all volumes. Left: all independent momentum 
components. Right: components near the diagonal (k, k, k, k), (k ± 1, k, k, k), k > 0. 



Figs U show a volume and lattice momentum cut dependence using the other representation of 
the one-loop dressing function Z^(aq) vs. (aq) 2 . This behaviour is similar in all loop contribu- 
tions that have been studied. In the left figure the different branches for off-diagonal momentum 
tuples are clearly seen, which do not possess a continuum limit a — > 0. Restricting to momentum 
values near the diagonal {k, k, k, k), {k ± 1, k, k, k), k > 0, a universal momentum dependence for 
larger volumes shows up. 

This universal curve can be compared with the known one-loop analytic result at L — > oo and 
a -► [@ 

Z {l \aq) = -0.24697 log (ag) 2 + 2.29368 . (5.1) 
The aim is to verify the constant 2.29368. We fit the dressing function near diagonal in the form 

Z {l \aq) = -0.24697 log(ag) 2 + C L + a{aq) 2 + c 2 (aqf (5.2) 

by assuming an universal anomalous dimension and parametrise additional lattice artefacts via 
coefficients c\ and c 2 . A typical fit at L = 16 is presented in the left Fig. ||. The results for different 




Cr 




Figure 3: Left: One-loop fitting at L=16, Cl=i6 
known infinite volume result Coo = 2.29368. 



2.3050(35). Right: C L vs. 1/L 3 together with the 



volumes together with the known infinite volume result is shown in the figure on the right. From 
here is no doubt that NSPT will reproduce the one-loop analytic result. 
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Finally we present in Fig. |] the perturbative dressing function Z(q, n max ) = Xm=cT (?) / ^ 
for near-diagonal lattice 4-tuples summed up to four loops for L < 16 and up to one loop for 
L = 16 using (3 = 6. 
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Figure 4: Perturbative dressing function 
^(9, "max) up to four loops (one loop) vs. q 2 
using ft = 6 at L = 8, 10, 12 (16). 



6. Summary 

In the present work we have applied NSPT to calcu- 
late the Landau gauge gluon propagator in higher- 
loop perturbation theory. Our results are in good 
agreement with expectations from standard LPT 
in one-loop. There is a good chance to extract 
higher loop finite contributions using known lead- 
ing and subleading anomalous dimensions. Our re- 
sults have to be confronted against non-perturbative 
Monto Carlo results and interpreted. The present 
work is in progress, and we will combine our efforts 
with those of the Parma group in order to study the 
perturbative gluon propagator at larger lattices and 
also the ghost propagator (not discussed here). 
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